Evidence of Vortex Jamming in Abrikosov Vortex Flux Flow Regime 
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We report on dynamics of non-local Abrikosov vortex flow in mesoscopic superconducting Nb 
channels. Magnetic field dependence of the non-local voltage induced by the flux flow shows that 
vortices form ordered vortex chains. Voltage asymmetry (rectification) with respect to the direction 
of vortex flow is evidence that vortex jamming strongly moderates vortex dynamics in mesoscopic 
geometries. The findings can be applied to superconducting devices exploiting vortex dynamics and 
vortex manipulation, including superconducting wires with engineered pinning centers. 
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The dynamic behavior of vortices in type-II supercon- 
ductors is important in applications [l| and provides an 
exemplary model system. Abrikosov vortices in type-II 
superconductors are set in motion when the Lorentz force 
due to local supercurrents exceeds the strength of vortex 
pinning forces, resulting in energy dissipation. There- 
fore, control of the dynamic behavior of vortices has a 
broad relevance for reducing losses in superconducting 
wires and microwave devices. Also, the ability to con- 
trol and manipulate single magnetic flux quanta could 
lead to novel superconducting, devices, such as current 
rectifiers 0-0] > logic elements [5, 6], or single photon de- 
tectors @- 

In this work we demonstrate vortex jamming effect 
through long-range manipulation of Abrikosov vortices 
in mesoscopic superconducting wires. Non-local vortex 
flow is initiated by applying a Lorentz force on parts 
of the vortex lattice, while observing vortex motion in 
areas where no external force is present. The nonlocal 
flow is due to the finite elasticity of the vortex lattice, 
i.e. finite tilt modulus C44 |8j, compression en and shear 
C66 moduli Enhanced rigidity of the vortex lat- 

tice (i.e. enhanced en and cqq) in mesoscopic super- 
conductors with reduced dimensionality (superconduct- 
ing strips) leads to the possibility of non-local manipula- 
tion of vortices on the order of several hundred vortex lat- 
tice spacings 0, [H, [Hj . Non-local vortex manipulation 
opens opportunities to sense perturbations of the vortex 
system far from the point where the external force acts 
on an individual vortex. This dramatically increases the 
sensitivity to dynamic properties of Abrikosov vortices in 
mesoscopic superconductors. 

Here we show that long-range manipulation of vortices 
in low- ft materials, such as Nb, can be exploited to study 
granular-like properties of vortex matter, including vor- 
tex jamming [14|,ll5|- Reduction or complete termination 
of vortex flow due to strong vortex- vortex interaction in 
constricted (funnel) geometries leads to a reduction of 
dissipation due to vortex movement. Similar to granular 
materials, the affinity to jam in funnel geometries is most 



pronounced at high particle densities, i.e. in our case at 
high magnetic fields. We examine the jamming effect in 
the case of a single constriction that is much larger than 
the size of a vortex. Vortices are channelled into a con- 
striction on the order of few vortex cores wide [15]. We 
show that the vortex mobility is asymmetric with respect 
to the direction of motion of the vortices, as observed by 
the nonlocal voltage in Nb bridges, and discussed theoret- 
ically by Vodolazov et al. [l2[ . The experimental evidence 
of granular-like behavior of vortex matter that we show 
in this paper should initiate wider application of theoret- 
ical concepts developed in the physics of granular mate- 
rials to the behavior of superconducting vortices. Based 
on these concepts, the engineering of pinning landscapes 
that suppress vortex motion via jamming will lead to 
lower dissipation at high magnetic fields in a wide range 
superconducting devices and wires. 

The 100-nm-thick Nb film was sputtered using a dc 
magnetron on a silicon wafer covered with 500 nm of sil- 
icon nitride. The critical temperature of the as-grown, 
high quality Nb film was 9.1 K with transition width of 
50mK and residual resistivity ratio RRRsimb. We used 
e-beam lithography and reactive ion etching (combina- 
tion of CF4 and SF6 plasma) to define structures simi- 
lar to those shown in the inset of Figure [U The 200-nm 
wide bridge had the normal state resistance of 8.65 /i(]cm. 
Considering that the slope dH c2 /dT near T c is 0.4 T/K 
and the fact that we are in the dirty limit, we obtain 
the basic Ginzburg-Landau superconducting parameters 
of the structure: the coherence length £(0) = 13.5 nm 
and the penetration depth A(0) = 145 nm, resulting in 
k= A(0)/1.63f(0)=6.6 PHI. 

The choice of material was motivated by the require- 
ment of strongly interacting vortex matter in order to 
reduce the compressibility of the lattice and enhance the 
vortex jamming effects. This can be accomplished by en- 
hancement of the elasticity moduli of the vortex lattice, 
which led us to use type-II superconductors with small 
values of ft = A/£ [17]. However, the intrinsic pinning in 
such thin films is strong due to the same reasons, and the 
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FIG. 1. (Color online) Local current- voltage characteristics of 
the Nb bridge at 2 K. Inset: atomic force microscopy image of 
the device used in the studies of nonlocal vortex movement - 
/- local current, Vi oc - local voltage, V n ii and V n i2 - nonlocal 
voltage leads. The distance between / and V n ii and V n i2 at the 
central vortex guide channel is 1.5 /im and 3 /im, respectively. 



vortex-pin interaction usually dominates over the vortex- 
vortex interaction. Therefore, the quality of the Nb thin 
film as well as precise definition of the mesoscopic struc- 
tures are both important. 

The Nb films were patterned into multi-terminal struc- 
tures, as shown in the inset of Fig. [TJ Direct current 
was applied through terminals /, and a local voltage was 
measured via terminals V\ oc (conventional four-terminal 
geometry), while the non-local voltage was measured via 
terminals V n i\ and V n i2 using a Keithley 2182 nanovolt- 
meter. The distance between the local leads to the V n \\ 
is 1.5 /im and to the V n \2 is 3.0 /im. The center supercon- 
ducting line with width w = 200 nm serves as a guiding 
channel for Abrikosov vortices. When the external mag- 
netic field is applied perpendicular to the plane of the 
film, vortices are nucleated in the superconducting pads 
located at either side of the sample, and depending on the 
direction of the current and the polarity of the vortices, 
the vortices can be moved either to the left or to the right 
side of the sample under the influence of Lorentz force. 
The coupling of the left and right superconducting Nb 
pads with the vortex guiding channel is through asym- 
metric constrictions. In case of granular-like behavior of 
the vortices, the funnel provides a well-defined guided 
concentrator for vortex movement to the right, facilitat- 
ing vortex jamming [14]. 

The geometry of the structure is such that the applied 
current is confined to the current lead / and exerts a 
direct force on the vortices that are only in that part 
of the structure. The applied current decays exponen- 
tially oc exp(— ttx/w) in the vortex guiding channel at 
a distance x away from the intersection between the lo- 



cal current lead and the vortex channel 0|. Therefore, 
no voltage is expected at both V n \\ and V n \2 when the 
structure is either in the normal or in superconducting 
state without vortex flow. 

Local current-voltage characteristics at 2 K in an ap- 
plied magnetic field show that we can isolate the vortex 
flux-flow regime in the range of currents between 100-200 
/iA (Fig. [J). The flux flow regime in this part of the vor- 
tex phase diagram is characterized by high vortex flow 
velocities [18|. The onset of the linear regime at ~ 100 
/iA can be interpreted as dynamic phase transition [l9| 
from plastic flow at I < 100 /iA to a coherent movement 
of the vortex lattice as a whole at I > 100 /iA. 

The magnetic field dependence of the nonlocal voltage 
is shown in Fig. [2^. The current of 150 /iA was ap- 
plied through the structure in order to initiate the flux 
flow through the superconducting vortex guiding chan- 
nel. The nonlocal voltage appears at nonzero magnetic 
field, reaches some finite value and then disappears at 
higher magnetic fields approaching H C 2 at this tempera- 
ture. The voltage signal shown in this figure is obtained 
by taking a difference in voltage signals for opposite di- 
rections of the applied current thus eliminating any ther- 
mal emf. The signal is symmetric about zero magnetic 
field. With increasing distance from the driving current 
the non-local voltage is reduced, but it is still finite even 
3 /im away. The non-local voltage decays with distance 
from the current lead where Lorentz force is present. If 
the vortex lattice was ideally rigid the non-local volt- 
age would not decay with distance. But finite rigidity of 
the vortex lattice causes local compression of the lattice 
and presence of disorder in the film causes finite pinning. 
Both of these contribute to loss of vortex flow along the 
vortex carrying stripe and decay of non-local voltage with 
distance. 

One stark difference from earlier reports on the nonlo- 
cal vortex motion is the appearance of highly symmetric 
peaks and kinks in the nonlocal voltage at specific val- 
ues of the magnetic field. The features are correlated 
between the first and second nonlocal probe, indicating 
highly correlated, long-range vortex flow along the chan- 
nel. The main features in the nonlocal voltage curves 
can be observed at 190, 435, 730, and 1070 mT. Weakly 
pronounced features, such as kinks, are present at inter- 
mediate fields, and they correlate between the two non- 
local voltage signals. Considering the mesoscopic size of 
the vortex channel the main peaks are to be associated 
with transitions between structural configuration of the 
moving vortex lattice in the channel. The narrow width 
of the channel that is on the order of the vortex- vortex 
spacing imposes its ID geometry on the vortex lattice. A 
triangular vortex lattice re-organizes into a vortex chain 
Since the driving current is strong, 
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such that the vortex lattice is moving with relatively 
high speed, the vortex lattice can be considered dynami- 
cally ordered [19] and its configuration can be calculated 
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FIG. 2. (Color online) (a) Nonlocal voltage as a function 
of applied magnetic field measured in the 200-nm-wide chan- 
nel at a distance of 1.5 /mi {Vnii - black) and 3 /mi (V n i2 - 
red) from the applied current. The curves are shifted ver- 
tically for clarity. Blue vertical dotted lines correspond to 
magnetic fields at which vortex chain transitions take place 
with schematic of vortex chain configurations at different field 
intervals. Arrows denote additional kinks in V n i(H) depen- 
dence; (b) magnetic field dependence of the vortex chain tran- 
sitions in the channel showing experimental values (points) 
and theoretical curve from Eq. (1). No fitting parameters are 
used. 



by applying equilibrium conditions. At fields below H2 
the vortices are arranged in a regular chain moving with 
steady velocity along the center axis of the mesoscopic 
channel. A Bean-Livingston barrier aligns the vortices 
in the center of the channel. When the magnetic field 
is increased, the intervortex spacing in the chain reduces 
gradually, until the intervortex distance becomes on the 
order of the half of strip width. At field H2 a structural 
transformation occurs - the vortex chain is suddenly split 
into two chains and the intervortex distance within each 
chain increases. This leads to a rapid increase of the total 
amount of vortices flowing through the channel leading to 
a peak in nonlocal voltage. Further increase in magnetic 
field preserves the two-chain configuration while grad- 



ually decreasing the intervortex spacing within each of 
the chains. The process of chain splitting repeats at H3 
and H4 with formation of three and four vortex chains, 
respectively with a corresponding peak in the nonlocal 
voltage. Ginzburg- Landau calculations [io| show that 
the magnetic fields at which the equilibrium reordering 
of the vortex lattice occurs are: 
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where A a 5 and A c are the in-plane and out-of-plane pene- 
tration depths, respectively, is the magnetic flux quan- 
tum, w is the width of the strip, and N=2,3,4... is the 
peak index. Since Nb film is in dirty superconducting 
limit, the anisotropy X a b/^c 1- The correspondence of 
the peaks to the theoretical model is shown in Fig. Et>. 
The slight discrepancy at higher order peaks might be 
due to the renormalization of the Bean-Livingston bar- 
rier at high fields, which would result in a slight change 
of the effective vortex channel width w. 

Additional more weakly pronounced kinks in the non- 
local voltage can be observed at magnetic fields that are 
between the H^r vortex transitions. In order to elucidate 
the nature of these additional features we go back to 
the origin of the nonlocal voltage. The nonlocal volt- 
ages V n ii and V n i2 shown in Figure [2^ are the aver- 
age absolute voltages induced by positive and negative 
driving currents 7 + and I~ in the current channel, i.e. 

V nl (I + ,H)-V nl (I-,H) 



V n l(H)-- 



The V n ii and V n i 2 signals 



are dominated by the vortex flux flow in the central chan- 
nel and any asymmetries of vortex motion up and down 
the channel are averaged out. In Fig. [3] we show the 
traces of the original V n i2(I + , H) and V n i2(I~ , H) from 
Fig. [2^. There are two distinct features in these curves. 
First, the signal is symmetric with respect to simultane- 
ous inversion of field and driving current, showing that 
there is an easy vortex flow direction that is independent 
on vortex polarity (vortices and antivortices flow easier 
in the same direction). Second, additional periodic kinks 
become more pronounced in the easy vortex flow direc- 
tion. 

The reason for the voltage asymmetry is due to the 
asymmetric geometry of the structure - the vortex guid- 
ing channel is coupled to the left and right endpoints 
through a different geometry. The superconducting pads 
that couple to the vortex guiding channel serve as vortex 
reservoirs that source vortices to the channel or drain the 
vortices from the channel. On the left end we have the 
funnel-shaped constriction and on the right end there is a 
square superconducting pad. The ability for the vortices 
to enter into the channel from these two pads is different; 
the funnel is designed to jam the vortices as they enter 
the ID channel. The situation corresponds to the closed 
end channel described in 12]. Jamming that occurs for 
vortices entering the channel at the funnel end, causes re- 
duction in vortex chain mobility in the channel and the 
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FIG. 3. (Color online) Non-local voltage V n i2(I,H) as a func- 
tion of the applied magnetic field for two opposite directions 
of the driving current (T=2.4 K). Vortices of opposite polar- 
ities exhibit the same rectification behavior when moving in 
the same direction. 

asymmetry of the V n \ with respect to the driving cur- 
rent. The same reduction will take place for vortices of 
opposite polarity (antivortices) subject to opposite driv- 
ing current, as observed in Fig. [3] The results are similar 
if one drives the vortices through the channel by apply- 
ing the Lorentz force with a supercurrent through the 
left-most pair of contacts shown in the inset of Fig. [3l 

Furthermore, the mobility of the vortex chains in the 
channel will be dependent on the matching of the vor- 
tex structure in the source/drain pads and the chan- 
nel. Matching of the vortex configurations in the chan- 
nel (vortex chains) to the vortex configuration in the 
pads (vortex lattice) would cause higher vortex mobil- 
ity (larger V n i ) due to easier entrance/exit of vortices, 
while non-matching configurations will cause reduction in 
V n i. In our experiments, when we are using the right pad 
(square) as a source of vortices going into the channel, 
we observe higher V n i . We also observe more pronounced 
periodic peaks in V n \ corresponding to the geometrical 
matching effect between the square vortex lattice in the 
pad and the vortex chain structure in the channel. Fur- 
ther theoretical work could, possibly correlate with these 
matching configurations. 

In conclusion, we have studied the dynamics of non- 
local Abrikosov vortex flow in mesoscopic superconduct- 
ing Nb channels. High rigidity of the vortex lattice in 
this material is responsible for quasi- ID long-range cor- 
related vortex movement. The magnetic field dependence 
of the non-local voltage induced by the flux flow shows 
that vortices form ordered vortex chains. Voltage asym- 
metry (rectifying effect) with respect to the direction of 



vortex flow is evidence that vortex jamming can signifi- 
cantly moderate vortex dynamics in mesoscopic geome- 
tries. The finding can be applied to future devices ex- 
ploiting vortex dynamics and vortex manipulation Q . 
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